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ABSTRACT 

The  amorphous  structure  of  four  CagoMgxZruo-x  alloys  (X  =  10,  15,  20,  and  25  at.%)  has  been 
modeled  using  a  set  of  experimental  neutron  and  X-ray  scattering  data,  peak  fit  analysis,  and 
Reverse  Monte  Carlo  (RMC)  simulation.  The  amorphous  structure  can  be  described  as  a  mixture 
of  Mg-  and  Zn-  centered  clusters,  with  Ca  dominating  in  the  first  coordination  shell  of  these 
clusters.  The  coordination  number  (CN)  of  10  (with  about  7  Ca  and  3  (Mg+Zn)  atoms)  is  the 
most  common  for  the  Zn-centered  clusters.  CN  =11  and  12  (with  about  7-8  Ca  and  4  (Mg+Zn) 
atoms)  are  the  most  common  for  the  Mg-centered  clusters.  Fife-fold  bond  configurations 
(pentagonal  pyramids)  dominate  (-60%)  in  the  first  coordination  shell  of  the  clusters,  suggesting 
densely  atomic  packing.  Bond  angle  distributions  suggest  the  near-equilateral  triangles  and 
pentagonal  pyramids  to  be  the  most  common  nearest  atom  configurations.  Two  Ca-Mg-Zn  alloys 
with  inferior  GFA  have  increased  amounts  of  (0,0,12,0),  (0,2,8, 1),  (0,2, 8, 2)  and  (0,4,4, 3) 
clusters,  which  suggests  that  these  clusters  may  facilitate  crystallization. 
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INTRODUCTION 


Recently  discovered  ternary  Ca-Mg-Zn,  Ca-Mg-Cu  and  Ca-Mg-Al  bulk  metallic  glasses  (BMGs) 
have  unique  properties,  which  distinguish  them  from  other  BMGs  [  1,2, 3, 4, 5].  They  are  based  on 
two  simple  elements,  Ca  and  Mg,  and  have  the  lowest  density  (in  the  range  of  1.6  -  2.4  g/cm  ) 
among  all  currently  known  BMGs  [6].  These  metallic  glasses  also  have  extremely  low  Young’s 
(-20-30  GPa)  and  shear  moduli  (-9-14  GPa),  comparable  with  the  moduli  of  human  bones  [7]. 
Some  of  these  ternary  alloys  have  extremely  good  glass  forming  ability  (GFA),  similar  to  GFA 
of  the  best  Pd-  and  Zr-  based  BMGs  [4,6,8]. 

Figure  1  shows  the  Ca-rich  comer  of  the  liquidus  projection  of  the  Ca-Mg-Zn  system  [9].  The 
compositions  of  the  ternary  Ca-Mg-Zn  ternary  amorphous  alloys  reported  in  [2,4]  are  marked  as 
solid  circles  and  their  critical  plate  thicknesses  are  shown  by  corresponding  numbers  (in  mm). 
The  presence  of  the  ternary  eutectic  reaction  provides  a  strong  liquidus  temperature  gradient  in 
the  selected  composition  area,  which  leads  to  high  sensitivity  of  GFA  with  alloy  composition. 

For  example,  the  maximum  plate  thickness  at  which  an  alloy  is  fully  amorphous  after  casting  in  a 
water-cooled  copper  mold,  can  be  increased  by  more  than  10  times,  from  0.5  mm  to  6-10  mm, 
within  only  5%  composition  range  [4],  Thermodynamic  analysis  of  the  onset  driving  force 
(ODF)  for  crystallization  of  different  phases  in  the  Ca-Mg-Zn  system  has  clearly  shown  [10]  that 
GFA  of  the  ternary  alloys  improves  with  a  decrease  in  the  ODF  of  competing  phases.  For 
example,  amorphization  of  Ca6oMgxZn40-x  alloys  competes  with  crystallization  of  CaZn  (at  X  < 
17.5at%)  and  CaMg2  (at  X  >  17.5  at.%)  phases  (Figure  2). 

Based  on  relative  atomic  sizes  [11]  and  the  efficient  solute-centered  cluster  packing  (ECP) 
model  of  amorphous  structure  [12]  the  Ca-Mg-Zn  ternary  glasses  are  represented  as  <10, 9> 
efficient  cluster  packing  structures,  where  Mg  solutes  are  surrounded  by  -10  nearest  neighbors 
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and  Zn  solutes  have  ~9  nearest  neighbors,  assuming  that  all  the  nearest  neighbors  are  solvent 
atoms  (i.e.  Ca).  Glass  compositions  predicted  from  the  ECP  model  range  from  Ca7oMgi5Zni5  to 
Ca53Mgi2Zn35  [13],  which  fit  within  the  trapezoid  shown  in  Figure  1. 

The  current  work  is  aimed  at  studying  the  amorphous  structure  of  a  series  of  ternary 
Ca6oMgxZn4o-x  (X  =  10,  15,  20,  25  at.%)  metallic  glasses  with  the  use  of  neutron  and  X-ray 
scattering  and  Reverse  Monte  Carlo  (RMC)  modeling  and  relating  this  structure  to  the  alloy 
compositions  and  GFA.  The  amorphous  structure  identified  in  this  work  is  compared  with  the 
amorphous  structure  predicted  by  the  ECP  model  [12], 

EXPERIMENTAL  PROCEDURES 
Sample  Preparation 

Four  amorphous  alloys,  Ca6oMgioZn3o,  Ca6oMgi5Zn25,  Ca6oMg2oZn2o,  and  Ca6oMg25Zni5,  were 
used  in  this  study.  Among  these,  the  Ca6oMgi0Zn30  and  Ca6oMg25Zni5  alloys  show  marginal 
GFA,  with  the  maximum  fully-amorphous  plate  thicknesses  achievable  during  water-cooled 
copper  mold  casting  of  0.5  mm  and  1  mm,  respectively  [2,4],  Two  other  alloys,  Ca6oMgi5Zn25 
and  Ca6oMg20Zn2o  are  better  bulk  glass  formers  and  their  maximum  fully  amorphous  plate 
thicknesses  were  identified  to  be  6  mm  and  4  mm,  respectively.  In  order  to  exclude  the  casting 
condition  effects,  the  amorphous  specimens  of  all  four  alloys  were  prepared  by  the  same  method 
of  melt  spinning  in  the  form  of  ~0.2  mm  thick  ribbons.  Prior  to  the  neutron  and  X-ray  diffraction 
analysis,  these  ribbons  were  ground  into  powder.  The  density  p0  of  the  amorphous  alloys  was 
measured  with  a  helium  pycnometer  AccuPyc  1330  VI. 03  and  the  values  (in  g/cm3  and 
atoms/A3)  are  given  in  Table  1.  The  following  formula  was  used  for  the  density  conversion: 

Po  [atoms/A3]  =  (NJY^CiA  0  p0  [g/cm3] ,  ( 1 ) 
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where  Na  is  the  Avogadro’s  number  and  c\  and  A  \  are  the  atomic  fraction  and  atomic  mass, 
respectively  of  the  element  i  (Ca,  Mg  or  Zn)  in  a  given  alloy. 

Neutron  and  X-ray  Scattering 

The  neutron  scattering  experiments  were  conducted  at  room  temperature  in  vacuum  using  a 
General  Materials  (GEM)  diffractometer  at  the  ISIS  high-intensity  pulsed  neutron  source 
(Rutherford  Appleton  Laboratory,  Didcot,  UK).  The  GEM  detector  has  eight  banks  to  collect  the 
data  within  the  Q  range  from  0.1  to  100  A'1,  where  Q  is  the  scattering  vector  of  length  4  ;r  sin  (9/A, 
for  a  neutron  of  wavelength  A  scattered  at  an  angle  20.  The  exposure  time  for  each  specimen 
was  about  210  min,  which  corresponded  to  the  accumulated  intensity  of  -5.7x1 06  scattered 
neutrons.  The  X-ray  scattering  experiments  were  conducted  at  room  temperature  in  argon 
atmosphere  on  Panalytical  Xpert-PRO  diffractometer  (ISIS  facility)  using  Ag  radiation.  The 
useful  Q  range  was  from  0.5  to  18  A'1. 

The  ISIS  developed  software,  Gudrun  and  GudrunX,  were  used  to  process  the  neutron  and  X-ray 
scattering  data  and  obtain  the  experimental  total  scattering  structure  factors  S(Q)  and  SX(Q),  as 
well  as  total  radial  distribution  functions  (RDFs)  G(r)  and  Gx(r).  The  superscript  V  indicates 
that  the  data  was  obtained  by  X-ray  scattering.  These  functions  are  related  to  the  partial  radial 
distribution  functions  (PRDFs),  gy(r),  as  follow  [14]: 

a  3  oo 

S  (0  =  X  cicAbj  J  Asij  (r)  -  l]sin  (Qr)dr  (2) 

U  ij= 1  0 

S’  (Q)  =  Y,  cfij  4^22^1]  r\ga  (. r )  -  l]sin(er)rfr  (3) 

Q  IJ-'  YCM®2“ 

i= 1 

G(r)  =  Y,  cicAbj  [§ij  <A)  - 1]  (4) 

i,j= 1 
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(5) 


1] 

,J"  2>,Z 

i= 1 

In  Equations  (2)  through  (5),  c;,  bi,  J](Q)  and  Z,  are,  respectively,  the  atomic  fraction,  coherent 
bound  neutron  scattering  length,  X-ray  scattering  factor  and  atomic  number  of  element  i  (Ca, 
Mg,  or  Zn),  p0  is  the  average  number  density  of  the  material  (in  atoms  per  A3),  and  r  is  a 
distance  from  an  arbitrary  point  in  the  amorphous  structure. 

Reverse  Monte  Carlo  Simulation 

The  three-dimensional  models  of  the  amorphous  structures  of  the  Ca-Mg-Zn  alloys,  which  agree 
with  the  experimental  data,  boundary  conditions  and  coordination  constraints,  were  generated 
with  the  use  of  the  Reverse  Monte  Carlo  (RMC)  simulation  technique  described  in  details  in 
[15], The  simulation  boxes  with  periodic  boundary  conditions  contained  from  1600  to  4000 
atoms.  The  input  data  used  in  this  simulation  were  the  experimental  neutron  and  X-ray  total 
structure  factors,  S(Q )  and  S"(Q),  and  radial  distribution  functions,  G(r).  The  main  coordination 
constraints  used  in  this  work  were  the  fixed  alloy  density  and  excluded  volume  (minimum 
interatomic  distances).  These  two  constraints  were  found  to  improve  the  separation  of  partials  in 
cases  where  the  separation  matrix  is  poorly  conditioned  [15,16],  Because  the  dominant  effect 
determining  the  amorphous  structure  is  the  atomic  packing,  included  information  on  the  atomic 
bond  distances  in  the  model  severely  limited  the  number  of  structures  that  were  consistent  with 
the  experimental  data.  The  bond  distances  were  estimated  from  the  available  data  for  CaMg2, 
CasZn,  Ca5Zn3,  CaZn,  and  MgZn2  intermetallics  [17].  Two  additional  coordination  constraints 
were  also  used  during  initial  stages  of  the  simulations,  Additional  assumptions  that  (a)  the 
average  nearest  bond  distances  between  six  pair  atoms  (Ca-Ca,  Ca-Mg,  Ca-Zn,  Mg-Mg,  Mg-Zn, 
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and  Zn-Zn)  are  the  same  in  four  studied  alloys  and  (b)  the  nearest  atom-pair  bond  distributions 
are  described  by  Gaussian  functions,  were  also  used  during  the  first  stages  of  simulation. 

The  RMC  computed  amorphous  structures  were  analyzed  with  the  use  of  special  programs, 
which  allow  calculation  of  PRDFs,  type  and  distribution  of  coordination  polyhedra  (atomic 
clusters),  neighbor  coordination  numbers,  triplet  angle  correlations,  etc. 

RESULTS 
Pair  Correlations 

The  eventual  matches  between  the  RMC  simulated  and  experimental  total  scattering  factors  and 
radial  distribution  functions  after  ~107  RMC  steps  are  demonstrated  in 
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Figure  3.  An  excellent  fit  indicates  that  the  RMC-simulated  amorphous  structure  describes  the 
experimental  spectra  very  well.  The  first  G(r)  peak  located  in  the  r-range  between  2.21  and  4.6 
A,  corresponds  to  positions  between  atoms  within  the  first  coordination  shell.  Six  PRDFs 
corresponding  to  the  simulated  amorphous  structure  are  shown  in 
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Figure  4.  A  single  gij(r)  peak  inside  the  first  coordination  shell  range  is  clearly  seen  in  Ca-Ca, 
Ca-Mg,  and  Ca-Zn  PRDFs,  and  the  bond  distances  ry  between  the  Ca-Ca,  Ca-Mg,  and  Ca-Zn 
pairs  corresponding  to  this  peak  are  identified  to  be  3.84  A,  3.48  A,  and  3.21  A,  respectively.  On 
the  other  hand,  the  Mg-Mg,  Mg-Zn  and,  especially,  Zn-Zn  bond  distances  within  the  first  shell 
(i.e.  inside  the  first  RDF  peak)  are  not  defined  by  a  single  g,j(r)  peak.  Instead,  several  g,j(r)  peak 
positions  are  common  (see 
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Figure  4).  For  example,  three  gMgMg(r)  peaks  are  seen  within  the  first  coordination  shell 
corresponding  to  the  characteristic  Mg-Mg  bond  distances  of  3.02  A,  3.54  A  and  3.80  A.  For  the 
Zn-Zn  PRDF,  four  peaks,  at  rMgMg  =  2.63  A,  3.05  A,  3.61  A  and  4.04  A,  are  present  within  the 
first  shell.  Similar  behavior  of  PRDFs  has  also  been  noticed  for  other  three  Ca-Mg-Zn 
amorphous  alloys. 

The  inter-atomic  bond  distances  calculated  from  the  position  of  the  first  peak  of  respective 
PRDfs  for  four  Ca-Mg-Zn  amorphous  alloys  are  given  in  Table  2.  These  calculated  ry  are  about 
2.3  -  4.1  %  smaller  than  the  respective  metallic  bond  distances,  rm,  and  about  6.5  -  9.8%  larger 
than  the  respective  covalent  bond  distances,  rc  (see  Table  3).  Using  the  lever  rule, 
fc  =  (rm  -  fij)/(rm  -  rc)  (7) 
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the  covalent  fraction  of  bonds,  fc,  was  determined  for  six  atomic  pairs  in  four  studied  amorphous 
alloys.  The  results  are  given  in  Table  4.  It  can  be  noticed  that  the  Ca-Zn  bonds  are  noticeably 
shorter  and  Zn-Zn  bonds  are  longer  in  marginal  glasses  (Ca6oMgioZn30  and  Ca6oMg25Zni5)  than 
in  good  bulk  glass  formers  (Ca6oMgi5Zn25  and  Ca6oMg20Zn20). 

Table  5  provides  characteristic  first-shell  bond  distances  between  Ca,  Mg,  and  Zn  atom  pairs  in 
several  binary  intermetallics.  The  crystal  structure  constraints  lead  to  discrete  values  of  the  bond 
distances  in  these  intermetallics.  It  can  be  seen  that  the  rCa-Mg  and  rZn-Mg  values  in  the  amorphous 
Ca-Mg-Zn  alloys  given  in  Table  2  are  smaller  than  the  respective  interatomic  distance  in  the 
crystalline  intermetallic  phases,  while  other  in  the  amorphous  alloys  are  close  to  the  minimum 
rij  values  in  the  crystalline  intermetallics.  This  observation  may  indicate  that  the  absence  of  long 
range  order  constraints  leads  to  shortening  the  interatomic  distances  between  the  nearest 
neighbor  atoms  in  the  amorphous  structure,  as  compared  to  the  binary  intermetallics. 

Weighted  bond  distances,  r"’e‘sh,ed ,  between  the  atom  pairs  within  the  first  peak  of  the  total  RDFs 

of  the  Ca-Mg-Zn  amorphous  alloys  are  given  in  Table  6.  They  are  calculated  using  the  following 
equation: 

=  )rgll(r)dr  }gj(r)dr  (8) 

'i  /  'i 

where  r  i  and  r2  are  the  minimum  and  maximum  cut-offs  for  the  first  shell.  When  the  values  in 
Table  6  are  compared  with  respective  values  in  Table  2,  it  can  be  found  that  r'p,ghted  are  always 

higher  than  the  first  peak  r,y  values.  The  difference  increases  with  a  decrease  in  the  size  of  pair 
atoms. 

The  clearly  identified  single  PRDF  peak  for  Ca-Ca,  Ca-Mg  and  Ca-Zn  atom  pairs  within  the  first 
coordination  shell  may  indicate  that  the  atoms  in  these  pairs  attract  each  other  and  are  touching 
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each  other.  On  the  other  hand,  Mg-Mg,  Mg-Zn  and,  especially,  Zn-Zn  pairs  are  less  attractive  to 
each  other.  As  a  result,  Mg-  and  Zn-  centered  atoms  prefer  to  have  Ca  atoms  as  the  nearest 
neighbors  in  the  first  shell  and  the  first-shell  Mg  and  Zn  atoms  just  fill  gaps  between  the  Ca 
atoms,  so  that  their  distances  from  the  center  atom  (Mg  or  Zn)  vary  depending  on  the  gap  values. 
This  observation  also  indicates  that  Zn  atoms  like  to  be  more  separated  than  other  atoms. 
Coordination  Polyhedra  and  Coordination  Numbers 

Atoms  in  metallic  glasses  tend  to  arrange  themselves  so  as  to  maximize  the  local  packing 
density[12,18  ,24]  and  the  finite  sizes  of  atoms  impose  geometrical  restrictions  on  the  possible 
local  atom  environments  [19].  The  local  environment  of  an  atom  i  is  generally  treated  in  terms  of 
the  geometry  of  the  shell  formed  by  its  nearest  neighbor  atoms  [23,20],  The  chemistry  (i.e.  the 
number  of  different  atom  species  and  their  relative  positions)  of  the  first  coordination  shell 
provide  additional  important  information  about  the  local  environment. 

Interatomic  bonds  in  metallic  glasses  are  generally  non-directional,  and  the  nearest  neighbors  for 
any  chosen  atom  i  are  often  defined  as  atoms,  which  are  in  contact  with  i  or  separated  from  it  by 
a  distance  smaller  than  a  certain  threshold  value.  The  latter  is  defined  globally  as  the  position  of 
the  first  minimum  of  RDF  of  the  whole  atomic  array  and,  therefore,  this  approach  is  not  sensitive 
to  the  details  of  the  first  shell  geometry  and  to  local  density  and  chemistry  fluctuations. 

Another,  more  advanced,  approach  is  to  define  the  nearest  neighbors  as  atoms  that  have  common 
faces  in  their  Voronoi  polyhedra  [18].  This  definition  of  the  first  coordination  shell  atoms  is 
straightforward  and  reflects  the  local  details  of  the  atomic  packing  topology.  In  this  approach,  a 
coordination  polyhedron  with  the  vertices  located  in  positions  of  the  first  coordination  shell 
atoms  and  with  the  edges  coinciding  with  the  interatomic  bonds  is  defined  for  any  chosen  atom  i 
[23],  Each  coordination  polyhedron,  which  is  also  called  an  /-centered  cluster  [12,21],  is 
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assigned  a  characteristic  Voronoi  signature  (v^w^.v^),  where  vm  is  the  number  of  the  vertices 
(first  coordination  shell  atoms)  with  m  entering  edges/bonds  [22]  (m  is  also  called  the  vertex 
coordination.)  The  /-centered  clusters  with  the  same  Voronoi  signature  are  considered  to  be 
topologically  equivalent  (even  if  they  are  not  identical)  because  they  can  be  transformed  into 
each  other  by  ‘elastic’  deformation  without  changing  the  number  of  all  vertices  and  connecting 
edges.  In  addition  to  assigning  the  topology,  the  Voronoi  signature  also  defines  the  total 
coordination  number  (CN)  of  the  cluster  as  CN  =  £vm.  The  same  Voronoi  signature  clusters  can 
however  be  chemically  different  if  they  consist  of  different  elements.  Therefore,  in  addition  to 
the  Voronoi  signature,  partial  coordination  numbers,  i.e.  the  number  of  atoms  of  different 
species,  should  also  be  known  for  a  more  complete  description  of  the  short  range  order  in  the 
amorphous  structure. 
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Figure  5  illustrates  types  and  fractions  of  Ca-,  Mg-  and  Zn-  centered  coordination  polyhedra 
found  in  the  amorphous  structures  of  four  Ca-Mg-Zn  alloys.  Although  many  types  of  the 
coordination  polyhedra  are  presented  in  the  amorphous  structures,  the  most  populated  clusters 
can  be  identified.  These  are  (0,1,10,2)  for  Ca-centered,  (0,2,8, 1)  for  Mg-centered  and  (0,3, 6,0), 
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Fraction 


(0,2,8, 1),  (0,3,6, 1)  and  (0,2, 8,0)  for  Zn-centered  clusters.  Three  examples  of  the  clusters 
extracted  from  the  amorphous  structures  are  given  in 


(a)  (b)  (c) 

Figure  6.  It  can  be  noticed  from 
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Figure  5  that  the  marginal  glass  forming  alloys  (with  10%  and  25%  Mg)  have  an  increased 
fraction  of  (0,0,12,0),  (0,2,8, 1),  (0,2, 8, 2)  and  (0,4, 4, 3)  clusters  in  comparison  with  the  good  glass 
forming  alloys  (with  15%  and  20%  Mg).  It  is  likely  that  these  clusters  facilitate  crystallization 
during  solidification  thus  reducing  glass  forming  ability. 
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Though  the  amorphous  structures  of  Ca6oMgxZn40-x  alloys  contain  many  types  of  the  first-shell 
clusters,  the  competing  crystal  phases,  CaMg2  and  CaZn,  [10]  contain  only  three  and  two  types 


of  the  clusters,  respectively.  For  example, 


(a) 


(b) 


(c) 


Figure  7  shows  three  principal  clusters  (one  is  the  Ca-centered  and  two  are  Mg-centered) 


presented  in  the  CaMg2  intermetallic  phase,  and  (a) 


(b) 
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Figure  8  shows  two  clusters  (Ca-  and  Zn-  centered)  presented  in  the  CaZn  intermetallic  phase. 
The  Ca-centered  clusters,  (0,0,12,4)  in  CaMg2  and  (0,1,10,6)  CaZn,  are  not  presented  in  the 
analyzed  amorphous  structures.  The  two  Mg-centered  clusters  in  CaMg2  are  (0,0,12,0)  and  the 
Zn-centered  cluster  in  CaZn  is  (0,3, 6,0).  Similar  Mg-  and  Zn-  centered  clusters  are  also 
populated  in  the  amorphous  structures  (see 
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Vofonoi  Polyhedra  Signature 


Figure  5). 

The  distributions  of  the  atomic  clusters  by  the  total  coordination  numbers  (CN)  in  the  first  shell 
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Fraction  Fraction 


Figure  7.  (a)  Ca-centered  (0,0,12,4)  and  (b,  c)  Mg-centered,  both  are  (0,0,12,0),  coordination 
polyhedra  in  the  CaMg2  crystal  structure. 


Figure  8.  (a)  Ca-centered  (0,1,10,6)  and  (b)  Zn-centered  (0,3, 6,0)  coordination  polyhedra  in  the 


CaZn  crystal  structure.  [] 
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Figure  9  for  the  Ca6oMgxZn40-x  amorphous  alloys.  There  is  no  notable  difference  in  the 
distributions  between  good  (X  =  15  and  20)  and  marginal  (X  =  10  and  25)  glass  formers.  The 
most  populated  coordination  numbers  are  CN  =  13  for  the  Ca-centered  clusters,  CN  =11  and  12 
for  Mg-centered  clusters  and  CN  =  10  for  Zn-centered  clusters.  The  CN  values  for  the  Zn-  and 
Mg-  centered  clusters  are  higher  than  the  CN  values  of  9  and  10,  expected  for  these  clusters  from 
the  densely-packed  solute-centered  cluster  (DPSCC)  model  [12].  This  discrepancy  is  most  likely 
due  to  the  fact  that  the  DPSCC  model  assumes  that  only  solvent  (Ca)  atoms  are  presented  in  the 
first  coordination  shell.  However,  if  some  Ca  atoms  are  replaced  by  the  smaller  Mg  and  Zn 
atoms,  larger  number  of  atoms  is  expected  in  the  first  shell. 
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Atom  Pair 


Figure  10  shows  the  average  number  of  Ca,  Mg,  and  Zn  atoms  in  the  first  shell  of  the  Ca-,  Mg-, 
and  Zn-  centered  clusters  in  the  Ca-Mg-Zn  amorphous  alloys.  It  can  be  seen  that  the  Ca  atoms 
prevail  in  the  first  shell  of  the  clusters  and  an  increase  in  the  concentration  of  Mg  from  10  to 
25%  almost  does  not  affect  the  number  of  Ca  atoms  in  all  three  kind  clusters.  The  respective  Ca- 
Ca,  Mg-Ca  and  Zn-Ca  partial  coordination  numbers  (PCN)  are  PCNca-ca  =  8. 1-8.5,  PCNMg-ca  = 
7.2-7. 9,  and  PCNZn-ca  =  6. 7-7.2.  At  the  same  time,  the  number  of  Mg  atoms  increases  (i.e. 
PNCca-Mg  increases  from  1.2  to  3.0  and  PNCMg-Mg  increases  from  0.9  to  3.7)  and  the  number  of 
Zn  atoms  decreases  (PNCca-zn  decreases  from  3.3  to  1.8  and  PNCMg-zn  decreases  from  2.3  to  0.6) 
in  the  Ca-  and  Mg-  centered  clusters  with  an  increase  in  the  concentration  of  Mg.  The  amount  of 
Zn  and  Mg  in  the  Zn-centered  glasses  is  mainly  dependent  on  the  type  of  glasses  than  on  the 
amount  of  Mg  in  the  alloys.  Namely,  marginal  glasses  (with  10%  and  25%  Mg)  have  an 
increased  amount  of  Zn  and  a  decreased  amount  of  Mg  in  Zn-centered  clusters  than  good  glass¬ 
forming  alloys  (with  15%  and  20%  Mg). 

The  neighbor  environment  of  a  specific  j  atom  in  the  first  coordination  shell  of  an  /-centered 
cluster  can  be  described  by  the  vertex  coordination  number  nij,  i.e.  the  number  of  the  j  atom 
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neighbors  located  in  the  first  shell  of  the  /-centered  cluster  [23,24]  or  the  number  of  edges 
(bonds)  entering  the  vertex  j. 


3-fold  4-fold  5-fold  6-fold  3-fold  4-fold  5-fold  6-fold 


Bond  Type 


Bond  Type 


3-fold  4-fold  5-fold  6-fold 
Bond  Type 

Figure  1 1  shows  the  fraction  of  3-,  4-,  5-  and  6-  fold  bond  vertices  in  the  Ca-,  Mg-,  and  Zn- 
centered  clusters.  In  all  three  cluster  types,  the  vertex  coordination  m  =  5  prevails.  Domination  of 
polyhedra  with  4-,  5-,  and  6-  bond  vertices  has  been  shown  [23,20]  to  be  indication  of  dense 
atomic  packing  of  the  amorphous  structure.  Moreover,  high  fraction  of  5-bond  vertices  can 
indicate  domination  of  pentagonal  bi-pyramids  as  structural  units  in  the  amorphous  structure.  It 
is  worth  to  note  that  two  pentagonal  bipyramids  with  a  common  5-fold  vertex  connected  by  ten 
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tetrahedrons  with  common  faces  form  a  coordination  icosahedron,  which  is  believed  to  be  most 
common  cluster  in  a  densely  packed  amorphous  structure  [25]. 

Bond  Angle  Distributions 

The  bond  angle  distributions  can  provide  additional  information  about  most  common  local 
structures  around  the  cluster  vertices,  because  characteristic  vertex  angles  are  very  sensitive  to 
the  vertex  type  (e.g.  number  and  length  distributions  of  vertex  bonds). 
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Figure  12  shows  Ca-Ca-Ca,  Ca-Mg-Ca  and  Ca-Zn-Ca  bond  angle  distributions  in  the  Ca-Mg-Zn 
amorphous  structures.  The  upper  limit  of  the  bond  lengths  has  been  set  to  4.6  A,  which 
corresponds  to  the  end  of  the  first  peak  of  the  total  RDFs.  A  characteristic  feature  for  all  four 
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alloys  is  the  first  peak  located  near  60°.  The  exact  60°  angle  position  for  the  Ca-Ca-Ca  triplets 
corresponds  to  the  close  packing  of  three  equal  hard  spheres  and  forming  equilateral  triangles. 

For  the  Mg-  and  Zn-  centered  triplets  (i.e.  Ca-Mg-Ca  and  Ca-Zn-Ca,  the  first  peak  is  located  at 
-64°  and  67-69°,  respectively.  These  shifts  of  the  first  peak  angle  position  to  the  higher  values 
can  be  explained  by  the  smaller  atomic  radii  of  Mg  and  Zn  atoms.  From  the  known  bond  angle 
values,  the  effective  (most  common)  bond  distance  ratios  are  estimated  to  be  r Ca-Mg/r ca-ca  =  0.94 
and  rca-Zn/i'Ca-Ca  =  0.89,  which  leads  to  the  effective  atomic  radius  ratios  of  (rMg/rca)ett=  0.88  and 
(j"zn/r c«)eff  =  0.78.  These  calculated  radius  ratios  are  higher  than  the  metallic  radius  ratios, 
(rMg/rca)m  =  0.81  and  (rZn/rca) m  =  0.68,  and  covalent  radius  ratios,  (rMg/rca) c  =  0.80  and 
(rzn/rca) c  =  0.69,  both  are  taken  from  Table  3.  The  observed  Ca-Mg-Ca  and  Ca-Zn-Ca  triplet 
configurations  can  be  explained  by  larger  shrinking  of  Ca-Ca  bonds  than  Ca-Mg  and  Ca-Zn 
bonds  (relative  to  the  respective  metallic  bond  distances),  which  is  in  agreement  with  the 
weighted  bond  distances  data  reported  in  Table  6. 

The  second  bond  angle  distribution  peak,  which  is  located  near  107°  for  Ca-Ca-Ca  triplets, 
corresponds  to  the  interior  angle  of  108°  of  a  regular  pentagon.  A  slightly  smaller  value  is 
probably  due  to  slight  distortion  of  the  pentagon,  when  one  of  the  five  Ca  atoms  in  the  pentagon 
vertices  is  replaced  by  Mg  and/or  Zn.  Due  to  shorter  Ca-Mg  (or  Ca-Zn)  bonds,  the  interior  angle 
at  the  Mg  (Zn)  vertex  of  the  pentagon  increases  at  the  expense  of  other  interior  angles,  the  total 
sum  of  which  in  the  pentagon  is  540°.  Indeed,  the  second  peak  for  Ca-Mg-Ca  and  Ca-Zn-Ca 
triplets  is  located  at  ~1 16°  and  -120°,  respectively.  Having  a  pentagon  consisting  of  4  Ca  atoms 
and  one  Mg/Zn  atom  and  using  these  values  of  the  interior  angle  at  the  Mg/Zn  vertex,  the  interior 
angles  at  the  other  four  Ca  vertices  are  calculated  to  be  106°/105°.  To  obtain  the  average  angle  of 
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107°  for  the  Ca-Ca-Ca  triplets,  it  should  be  assumed  that  one  Mg/Zn  atom  replaces  every 
tenth/fifteenth  Ca  atom  in  the  pentagon  configurations  presented  in  the  amorphous  structure. 

The  two  characteristic  peaks  in  the  bond  angle  distributions  that  correspond  to  the  angles  at  the 
vertices  of  equilateral  triangles  and  pentagons  support  our  earlier  conclusion  that  tetragons  and 
pentagonal  bipiramids  are  the  most  common  atom  configurations  in  the  Ca6oMgxZn40-x 
amorphous  alloys. 

CONCLUSIONS 

1 .  Using  combination  of  neutron  and  X-ray  scattering  and  Reverse  Monte  Carlo  modeling 
allowed  us  to  model  the  amorphous  structure  of  Ca6oMgxZn40-x  alloys  and  extract  partial 
distribution  functions  (PDF),  which  describe  well  experimental  neutron  and  X-ray  structure 
factors  and  total  correlation  functions  of  these  alloys. 

2.  Amorphous  structure  of  Ca-Mg-Zn  BMGs  can  be  described  as  a  mixture  of  Mg-  and  Zn- 
centered  clusters,  with  Ca  dominating  in  the  first  coordination  shell  of  these  clusters. 

3.  CN  =10  [~7  Ca  +  3  (Mg+Zn)]  is  most  common  for  Zn-centered  clusters. 

4.  CN  =11  and  12  [~7-8  Ca  +  4  (Mg+Zn)]  are  most  common  for  Mg-centered  clusters. 

5.  5-fold  bond  configurations  (pentagonal  pyramids)  dominate  (-60%)  in  the  first  coordination 
shell  of  the  clusters,  suggesting  densely  atomic  packing. 

6.  Bond  angle  distributions  suggest  the  near-equilateral  triangles  and  pentagonal  bi-pyramids  to 
be  the  most  common  atom  configurations. 

7.  Two  Ca-Mg-Zn  alloys  with  inferior  GFA  have  increased  amounts  of  (0,0,12,0),  (0,2,8, 1), 
(0,2, 8, 2)  and  (0,4, 4, 3)  clusters,  which  suggests  that  these  clusters  may  facilitate  crystallization. 
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TABLES 


Table  1.  Density  (in  g/cm3  and  atoms  per  A3)  of  Ca-Mg-Zn  amorphous  alloys  produced  in  this 
work. 


Alloy 

Density 

g/cm3 

A^ 

Ca6oMgioZn3o 

2.4481  ±0.0058 

0.031985 

Ca6oMgi5Zn25 

2.2890  ±  0.0073 

0.031301 

Ca6oMg2oZn2o 

2.1499  ±0.0024 

0.030837 

Ca6oMg25Zni5 

2.0043  ±  0.0044 

0.030227 

Table  2.  The  interatomic  bond  distances  (r;/,  in  A),  which  correspond  to  the  position  of  the  first 
peak  in  respective  partial  PRDFs  of  the  Ca6oMgxZn4o-x  amorphous  alloys.  The  bond  distance 
values  averaged  for  four  alloys,  (r;y)ave r,  as  well  as  the  deviations  of  these  values  from  metallic, 
rm,  and  covalent,  rc,  bond  distances  (see  Table  3),  are  also  given  here. 


rip  A 

V Ca-Ca 

f Ca-Mg 

F Ca-Zn 

7* Mg-Mg 

F Mg-Zn 

^Zn-Zn 

Ca6oMgioZn3o 

3.82 

3.44 

3.14 

3.02 

3.02 

2.66 

Ca6oMgi5Zn25 

3.80 

3.48 

3.19 

3.06 

2.78 

2.56 

Ca6oMg2oZn2o 

3.84 

3.48 

3.21 

3.02 

2.84 

2.56 

Ca6oMg25Zni5 

3.81 

3.52 

3.16 

3.09 

2.85 

2.63 

ir  ij)  aver 

3.82 

3.48 

3.18 

3.05 

2.87 

2.60 

\{){)%(raver/rm-\) 

-3.1 

-2.5 

-4.1 

-4.8 

-2.3 

-2.9 

1 00%(raver/rc- 1 ) 

8.5 

9.8 

6.5 

8.1 

9.2 

6.7 
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Table  3.  Metallic,  rm,  [26]  and  covalent,  rc,  [27]  bond  distances  (in  A)  between  Ca,  Mg  and  Zn 
atom  pairs. 


Ca-Ca 

Ca-Mg 

Ca-Zn 

Mg-Mg 

Mg-Zn 

Zn-Zn 

T m  (^) 

3.94 

3.57 

3.31 

3.20 

2.94 

2.68 

rc  (A) 

3.52 

3.17 

2.98 

2.82 

2.63 

2.44 

Table  4.  The  covalent  fraction,  A,  of  Ca-Ca,  Ca-Mg,  Ca-Zn,  Mg-Mg,  Mg-Zn,  and  Zn-Zn  bonds 


in  four  Ca-Mg-Zn  amorphous  alloys. 


fc 

Ca-Ca 

Ca-Mg 

Ca-Zn 

Mg-Mg 

Mg-Zn 

Zn-Zn 

Ca6oMgioZn3o 

0.29 

0.33 

0.52 

0.47 

0.00 

0.08 

Ca6oMgi5Zn25 

0.33 

0.23 

0.36 

0.37 

0.52 

0.50 

Ca6oMg2oZn2o 

0.24 

0.23 

0.30 

0.47 

0.32 

0.50 

Ca6oMg25Zni5 

0.31 

0.13 

0.45 

0.29 

0.29 

0.21 

(/c)aver 

0.29 

0.23 

0.41 

0.40 

0.28 

0.32 

Table  5.  Characteristic  bond  distances  in  several  binarv  intermetallic  nhases. 


MgZn2 

3.17 

3.04 

2.53 

2.61 

2.64 

Ca5Zn3 

3.61 

3.14 

2.69 

3.76 

3.29 

4.07 

3.34 

4.19 

Minimum 

3.61 

3.62 

3.14 

3.05 

3.03 

2.53 

Maximum 

4.32 

3.62 

4.17 

3.18 

3.04 

4.15 

Table  6.  Weighted  bond  distances,  r.eighted  (in  A),  between  the  atom  pairs  within  the  first  peak  of 

the  total  RDFs  (the  first  coordination  shell)  of  the  Ca6oMgxZn40-x  amorphous  alloys.  The  bond 
distance  values  averaged  for  four  alloys,  as  well  as  the  difference  between  the  calculated 
weighted  and  the  first  peak  bond  distances  (taken  from  Table  2)  are  also  given  in  two  last  rows. 


weighted  % 

rij  > A 

Ca-Ca 

Ca-Mg 

Ca-Zn 

Mg-Mg 

Mg-Zn 

Zn-Zn 

Ca6oMgioZn3o 

3.83 

3.55 

3.32 

3.35 

3.20 

3.21 

Ca6oMgi5Zn25 

3.85 

3.57 

3.29 

3.40 

3.33 

3.30 

Ca6oMg2oZn2o 

3.85 

3.53 

3.35 

3.44 

3.34 

3.36 

Ca6oMg25Zni5 

3.85 

3.52 

3.43 

3.43 

3.19 

3.38 

Righted  ayerage 

3.84 

3.54 

3.35 

3.40 

3.26 

3.31 

weighted  % 

rn  -  rip  A 

-0.10 

-0.03 

0.04 

0.20 

0.32 

0.63 
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FIGURES 


Mg 


Figure  1.  Liquidus  projection  of  Ca-Mg-Zn  ternary  system.  The  projection  shows  dashed 
isothermal  lines  and  corresponding  liquidus  temperatures  (in  °C)  and  thick  solid  lines  bounding 
Ca,  CaMg2,  CaZn,  and  CasZn  phase  fields.  Ternary  eutectic  and  peritectic  compositions  are 
marked  by  an  open  triangle  and  an  open  circle,  respectively.  The  compositions  of  amorphous 
alloys  and  their  maximum  plate  thicknesses  (in  mm)  are  shown  as  solid  circles  and  pertinent 
numbers.  A  trapezoid  with  a  dashed  boundary  represents  a  composition  range  for  the  ternary  Ca- 
Mg-Zn  BMGs,  in  accord  to  [4], 
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Onset  driving  forces,  ODF/RT 


2.5 


Figure  2.  Calculated  onset  driving  forces  of  various  crystalline  phases  for  Ca6oMgxZn4o-x  alloys 
versus  Mg  content  at  390  K.  The  circles  represent  the  amorphous  alloy  compositions  with 
indicated  critical  thicknesses.  [10] 
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Figure  3.  Experimental  and  RMC  simulated  X-ray  and  neutron  total  scattering  structure  factors 
S(Q)  and  neutron  total  radial  distribution  function  G(r)  for  a  Ca6oMg20Zn2o  amorphous  alloy. 
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Figure  4.  Partial  radial  distribution  functions  g(r)  for  Ca-Ca,  Ca-Mg,  Ca-Zn,  Mg-Mg,  Mg-Zn  and 
Zn-Zn  atomic  pairs  in  a  Ca6oMg20Zn2o  amorphous  alloy. 
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Figure  5.  Types  of  Ca-,  Mg-,  and  Zn-  centered  clusters  and  their  fractions  in  four  Ca6oMgxZn40-x 
amorphous  alloys. 
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Figure  6.  Configurations  of  three  typical  (a)  Ca-centered  (0,1,10,2),  (b)  Mg-centered  (0,2,8, 1), 
and  (c)  Zn-centered  (0,3,6, 1)  clusters  in  the  Ca6oMgxZn4o-x  amorphous  alloys. 


Figure  7.  (a)  Ca-centered  (0,0,12,4)  and  (b,  c)  Mg-centered,  both  are  (0,0,12,0),  coordination 
polyhedra  in  the  CaMg2  crystal  structure. 


Figure  8.  (a)  Ca-centered  (0,1,10,6)  and  (b)  Zn-centered  (0,3, 6,0)  coordination  polyhedra  in  the 
CaZn  crystal  structure.  [] 
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Fraction  Fraction 


Coordination  Number 


Coordination  Number 


Figure  9.  Distributions  of  Ca-centered,  Mg-cantered,  and  Zn-centered  clusters  by  their  total 
coordination  numbers  in  four  Ca6oMgxZn40-x  amorphous  alloys. 
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Figure  10.  Partial  coordination  numbers  in  Ca6oMgxZn4o-x  amoprhous  alloys  as  functions  of  Mg 
content. 
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Figure  11.  Fraction  of  vertices  with  3,  4,  5,  and  6  neighbor  atoms  (bonds)  in  Ca-,  Mg-,  and  Zn- 
centered  clusters  in  four  Ca6oMgxZn4o-x  amorphous  alloys. 
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Figure  12.  Ca-Ca-Ca,  Ca-Mg-Ca  and  Ca-Zn-Ca  bond  angle  distributions  in  Ca6oMgxZn4o-x 
amorphous  alloys.  A  pentagonal  pyramid  with  two  characteristic  peak  bond  angles  is  shown 
the  insert  and  represents  the  most  common  configuration  in  the  first  coordination  shell. 
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